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MEDITATION ON THE IDEA OF PONCELET’S THEOREM. 
[Philosophical Magazine, xx. (1860), pp. 8307—316.] 


HITHERTO Poncelet’s theorem has been regarded as a method sui generis 
and complete in itself; but in truth it is but the first germ or rudiment of a 
vast and prolific algebraical theory ; and not only so, but the principle which 
it contains admits of applications of the utmost value in various dynamical 
and analytical questions, which it is surprising should have been allowed to 
lay so long dormant. For the present, however, I mean to confine myself to 
a very brief indication of one direction in which the theorem admits of being 
generalized. And first I will make a remark upon so simple a matter as the 
extraction of the square root, which seems to have escaped observation, 
and at all events is so far from being generally known, that two of the 
highest authorities for mathematical erudition in this country whom I have 
consulted on the subject provisionally accept it as new. 


Let r be an approximate value of yN ;-then by that mode of application 
of Newton’s method of approximation to the equation a° = NV which is equi- 
valent to the use of continued fractions, we may easily establish the following 
theorem, namely, that 


CAN r+3rN m+ 6r°N + N* r+ 10° N + 5rN* 


2r ' +N’  4e+4eN ? Brel + NP?“ 


will be successive approximations to VN, whose limits of error can be 


* 


* In other words, if r be the first approximation to ,/N, the ith approximation will be 


(r+./N)§+(r— JN) 
(r+JN)— (r= Jyh’ 


so that the relative error becomes 

DDE ey 

(r+ /N)f—(r- JN) 
in which form the theorem is self-subsistent, and needs no proof. But the fact remains interest- 
ing, that the application of Newton’s method of approximation to the equation z?=N will be 


found to lead to the form above written at the ith step of the process conducted after the con- 
tinued-fraction fashion, 
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assigned when a limit to the error of the first. approximation r is given. The 
coefficients of the qth approximation, it will be observed, are for the numerator 
the alternate binomial coefficients 


—1 —-lq-2q-3 
uF 1L ? gi 1-4, be; 


and for the denominator the intermediate ones, 


r+ 3r N 
Br+ NN’ 
special case of a formula for any root of N given in the books; and to 
Mr De Morgan I am indebted for a hint which has led me to notice that all 
these forms may be deduced from the Newtonian method of approximation *. 


Mr Cayley has reminded me that the third approximation, is a 


If we call the ith approximation ¢ (7, r), we shall find that the functional 
equation ¢ {j, (i, r) =¢ (ij, r) will be satisfied; which is not so mere a 
truism as might at first sight be supposed, as any one may satisfy himself by 
studying the analogous theory for cubic or higher roots, a part of the subject 
to which I may hereafter return. 


Now as to the limits of accuracy afforded by the successive approxi- 
mations. Let e be a known limit to the relative error of the first approxi- 
VN—-r 

VN 
I take separately the cases where r is too great and r is too small. 


2 
mation r, by which I mean that ( ) <e. For greater simplicity, 


l. Let /N<r<(1+e) 7X; then the errors will be throughout in 
excess; and we may assign as a limit of error to the 7th approximation a 


quantity, say e;, which is a known function of e, namely, , which 


2 
(2e741)'-1 
s ' ; & 
it may be noticed is less than zm 

2. Let /N>r>(1—7)/N; then the errors will be alternately in defect 
and excess, and to the 7th approximation we may assign a limit of error ni, 
where n= -————__..—_-_.. f. 

™= E S TT e y A 

* The expansion (after Newton) of ,/N introduces the binomial coefficients—a curious fact! 
What are the analogous integers which the continued-fraction process applied to ¥/N will produce? 
+ If we write 
e;=0 (e, i) and n;=9 (n, i), 
then if i be any odd number, 
0 {0 (e, i), J} =O (6, ti), 
S {9 (m, i) j} => (m Ü); 
and if i be any even number, 
> {0 (e, i), J} =ð (e, Ùj) 
0 (> (n, i), J} =% (m, ij). 
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We may now apply these results to Poncelet’s linear approximate repre- 
sentation of y(a +bs + cx’), Suppose f+ ge is the first approximation, as 
found by Poncelet’s method, with a maximum relative error e, then 


(f+ gay + (a+ be + ca?) 
2 (f +ga) 
will be a much closer approximation, with a relative error never exceeding 


eet $ ; ; > j 
gpg in excess, nor so in defect. So a still nearer approximation will be 


(f+gzy +3 (f+gzx) (a+ ba + ca’) 
3(f+gaye+atbe + cx 


, with a relative error never exceeding 


Isiiri in excess, nor ieee in defect, and so on. The marvellous 
facility which these formule afford for the calculation of elliptic and ultra- 
elliptic functions, and not merely for their computation as by a method of 
quadratures, but (which is of far greater importance) their quasi-repre- 
sentation under circular and logarithmic forms, with assignable limits of 
proportional error, will be illustrated in a future communication. As regards 
the idea of substituting rational for irrational functions, I have only to-day 


learned from Mr Cayley that I am anticipated in this by Mr Merrifield*, 


Or more simply, if the error in excess be treated as positive, and in defect as negative, and 5 be 
the first and ô; the ith limit of error, we shall have 


 Wehack oa 
i— (2+8) - ôi’ 
and calling 6;=6 (i, ô), 
8 {j, 0 (i, 6)} =@ (ij, ô). 
Thus, then, if we call hae yx, yx will correspond to the (2%)th order of approximation, 
and the absolute value of the error will be less than 


Oe as 
(2 + 6)** — 5% 
By way of example, suppose we take 6 as our first approximation to ,/31, then 
cis i i if i 
and if we make patte, we shall have 
6 : 4/81 ::1+w:1, 
where gei 


2816.1’ 
which serves to exemplify the prodigious rapidity of the approximation in this method of extract- 
ing the square roots of numbers, 

* I quite concur with Mr Merrifield, and in fact before being made acquainted with the 
existence of his paper, had emitted the same opinion (among others to Dr Borchardt of Berlin), 
that the substitutive method, consisting in the employment of rational functions in place of the 
radical, affords by far the most expeditious means for the calculation of elliptic functions of 
all orders, especially the third, and supersedes the necessity for the construction of special 
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in a paper very recently read before the Royal Society, but not yet printed 
in the Transactions*. 


auxiliary tables. I believe, however, that my substitutions, founded on Poncelet’s views, are in 
general the best that can be employed for the purpose. In addition to other advantages they 
possess this, which deserves notice—that as we know œ priori a superior limit to the proportional 
error, the arithmetical values of the integrals to which they are applied may be brought out 
correct to any required place of decimals, without its being necessary to calculate and compare a 
superior and inferior limit to the integral, either one of these being sufficient in my method to 
indicate its own reliable degree of precision. 

* In general it is obvious, if px between the limits a and b retain always the same sign, 
and yx within these limits be sometimes greater and sometimes less than ¢x, but the difference 


a 
between them be always less than eg, then | “da wx will differ from | dx px by considerably less 
A b b 
a 
than e | dx@px. Paradoxical, however, as it may at first sight appear, there are extreme cases 
b t a 


where this difference tends to a ratio of equality with e "dx gx. The complete elliptic function 
b 


of the first order may be made to furnish an example of this. Let 


pare etay shoal napa bln aA tet) 
Ni-a?) (L-c?a?)} (1-a) (1 - a?) 
(so that b2=1~), and let 


_f¥(L-2*)+gbx , 
r= 1 — a?) J 2) } 


if we make 
2 TE. akl 
STEIF’ ea EAT TE K 


it follows from Poncelet’s theorem, that for all values of æ intermediate between 0 and 1, yx will 
differ from øx by less than epz. 


Now it will easily be found by ordinary integration that 
1 y Us Aa: Pe 


1 
Hence | dx gx must be always less than 
0 


ee ASTIE VONE X: ACNEE Bee eee eA 
fade vie (hee cL 

: 1 l+e 1 b 

t re Se S a 
that is, <5 log ;—, + ; tan me? 


when c becomes indefinitely near to unity; that is, when b becomes indefinitely small, this 
e Fai 
b ; 


approximate value for the integral in such case is log 5 so that the superior limiting value of 


approaches indefinitely near to log But we know, by a theorem of Legendre, that the 


1 
J dx x, found by the application of Poncelet’s method, approaches in this instance indefinitely 
0 


near to the value itself. The explanation of this is easy. As c approximates to unity, the only 
important values of x in the integral 


} dx uy 

o/{(1 = 22) (= ey 
are those which lie in the immediate vicinity of 1; and for all such values the relative error is at 
` a negative maximum. 
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The method, however, of Mr Merrifield in working out this conception is, 
I believe, entirely different from that here indicated: how the many mathe- 
maticians of a practical stamp, English and foreign, who have worked with 


It is not a little remarkable that so rude an application of Poncelet’s method should serve to 
indicate almost with the force of rigorous demonstration the approximate formula 


F (c)=log i + constant, 


when c approaches indefinitely near to unity, the constant left undetermined being known to be 


less than log 245. 


Nay, the demonstration may be made absolutely rigid if we set about to find an inferior 


limit. To this end make 
1 


POE CFL a8) + obo} J-a) 
we shall find without difficulty F 
1 io l+y-by+b Je 
I; aat jogir PeT A where y=4/(1+b°?), 
and consequently we shall obtain as an inferior limit to F (c) the expression 


1 l+y-by+b 
y °y-l+by-b’ 


which approaches indefinitely near to log > as c approaches indefinitely near to = It is thus 
seen that Legendre’s F (c), when c is itt near to 1, lies between log $ = and log $ +a 
the arithmetical mean between these limits is log 2 bt i , that is, log i+ 1:4785, differing by only 
-0923 from the true value log 7 + log 4. Of course, when the form of F (c) in the case supposed is 


known, namely, log; +0, there is no difficulty in determining C (as may be seen in Verhulst’s 


Traité des Fonctions Elliptiques) ; but the process above given of throwing the general value of F (c) 
between limits, is, I believe, by far the easiest and most natural method of obtaining this form. 
The limits themselves, it should be noticed, have virtually been found by the method, simple to 
naïveté, of writing ,/(1—c*x*)=,/(p?+q"), where p=,/(1-—2*) and q=ba, and then substituting 
s ki - as an inferior, and S 

* It p+a p*+q° 
Closer and calculable limits ad libitum to the integral may be arrived at by substituting for 


1 
Jt one or the other of the two following rational functions of p, q, according as we wish 


as a superior limit in the quantity to be integrated. 


to obtain an inferior or superior limit to the integral, namely, 
ipt+q+/(m+9")}'-{p+q-V(P+q)}i 1 


(P+ dt/(P +I (pt a— (P+ Vg) 
Pe sete Ae preat 1 

tP +I +HNAIP* +g) pta- Np +a’) Vg)’ 
in which formulæ the greater i is taken the closer will be the approximation. Iam not aware 
that any of these limits to F (e) (even the simplest of which, namely, those given above, may have 


some value for computational purposes, and have fallen thus very incidentally in my way) have 
ever before been noticed. 


It is not unworthy of notice that the second superior limit to a = —, namely, 
+ 
P+pa+¢? ih 


(p+a(P'+9)’ is an arithmetic mean between the first superior and first inferior limits, and 
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Poncelet’s method during the last quarter of a century, should have managed 
to overlook so obvious and important an extension of the principle and its 
applications, I find hard to realize; and my wonder is even greater that 
I should not have been anticipated twenty years ago, than that I should 
have been anticipated so recently. But the algebraical theory to which this 
extension points the way is replete with interest of a far higher order than 
its applications to practice; for plainly the derived approximate fractions, 
however sufficient for the purposes of computation, are not, nor ever can be 
the best and closest of their respective kinds*. To fix the ideas, let us 


consequently our second superior limit to the integral when b is indefinitely small becomes 
T 
4 
limit ; and as this approximation will evidently not stop at the second step of the process, we 
may safely infer that the integral derived from either formula when i= (for all values of b, 
whether finite or indefinitely small), not merely bears to F (c) a ratio differing infinitely little 
from that of equality, but is absolutely equal to, and may for all analytical purposes be employed 
to represent F (c). 

I have been at the trouble of calculating the inferior limit afforded by the second approxima- 


log - + -, which brings the constant much nearer to its true value than did the use of the first 


tion, and find that for b indefinitely small it is log : = -that is, log ++ 1:2977; the superior 


T 
3,/3 , 
limit has been shown to be log } + 14785, the mean is therefore log ; +1°3881, differing by only 
‘0018 from the true value! As the constant continues for all values of i to be a multiple of r, 
the ith approximations à supra and ab infra, which are always effectible, will give (on making 
i= ) two new expansions for 7, one infinitesimally in excess, the other infinitesimally in defect 
of its true value expressed as a multiple of log 2, which it might well repay the trouble of some 
young analyst to develope. 

* That the fractional forms derived from the linear substitutive form are not the best of their 
respective kinds, appears immediately, so far as the derivatives of the odd order (subsequent to 
the first) are concerned, from the consideration that the limits of error in excess and in defect 
will be actually attained for values of x lying within the prescribed limits; but these errors, 
e; and n; (when e=y, which is true by hypothesis), are never equal, the former (the extreme error 
in defect) being always the greater of the two; but if any such derivative were the best of its 
kind, the absolute values of the extreme errors of excess and defect ought to be equal to each 
other. But more generally, if possible, let the ith derivative to L (x) (where L (x) represents the 
radical linear approximant fx+b to ,/(a+bx+cx*), say Q(x)), namely, 

ee ee i) 
(La + Qx)’ — (Lx - Qr)? ; 


2 (La - Qx)? 
(La + Qa)*— (La — Qx)i 
maximum value of this must be equal (to the sign prés) to the value which it has when we give to 
x either of its extreme connecting values. Now obviously the above is a maximum only when 
Lux + Qx Lx 


- is a minimum, and therefore when — 
Lx- Qx z Qx 


x, say m, which makes this a maximum, gives to ma 1 the same value with the opposite sign 


be supposed the best of its kind: then the relative error is , and the 


is a maximum ; but by hypothesis, the value of 


to that which it would have in writing for x either of its limiting values, say k or k’. 


Thus we have two equations for determining Ok > Om’ namely, 
Lm Lk 
fe dain n 
Qm 3 Qk’ 
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confine ourselves to the second Ponceletic approximation to /(a+ bx + ca*), 
A+ wae + va 
1+ qa 
mined. The problem to be solved is the following. 
Let A+ prt vaer=V, 
(1 + qx) (a+ ba + ca*) = U; 
it is required to assign the four constants AÀ, p, v, q, so that the maximum 


namely, that which has the form : where 2, p, v are to be deter- 


value of u -1) for values of æ intermediate between a and b shall be the 


U 
least possible. Some little way, but only a little way, into the solution of 
this problem we can look in advance. In the first place, if we seek for the 


maximum values of (7 — 1) , we obtain the rational equation 


U 
dV dU 
Mat be + ca) (UT — v<.)- 0, 


‘which will easily be seen to be a cubic (not a biquadratic) equation in æ. 


U 
T, £a, £, the law of equality explained in my preceding paper would seem to 
show* that we must be able to satisfy the following equations, 


(por, = (per) = (Pa, = (pay? = (PbY, 
which amount to four independent equations, the precise number of constants 
A, 4, v, q to be determined. So in like manner the ith rational approximant 
will contain 2i disposable constants; the differentiation of the quantity 


Call (y — 1)=¢ (2); then the three roots of this equation being named 


analogous to 7 will give rise to an equation of the (2i—1)th degree; and 


wm EN-E a- (Gh) 


Im- 2k , which is incon- 


Qm Qk 


sistent with the first; so if i=3, we should obtain (3) = (a) , and therefore, on account 
L(m) 


of the first equation, Q E 1; and so in like manner for any value of i, we should derive one or 


Thus, suppose i=2, we should obtain from the second equation 


L a 
more numerical values for = which is absurd, since this quantity is a function of k, k’, the two 


connecting values of x. 
* Is it not, however, somewhat uncertain whether the equalities 


(G2 )?= (pæ) = ( gary)? 
must all, in all cases (that is to say, for all given values of the limits) subsist? since the law 
of equality will not apply to such values of x as lie without the prescribed limits, and non constat 
à priori that the roots of the cubic do all lie within these limits. The subject at the very 
threshold is beset with doubts and difficulties of a peculiar kind, which we can hardly hope to 
overcome without calling in geometrical imagination to our aid. 
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there will be 2c — 1+ 2, that is, 2i+ 1 functions of these 2i quantities to be 
equated, which furnish precisely the required number of equations to make 
the problem definite. It is, however, apparent that in solving these equa- 
tions we shall find a multiplicity of systems, by which I mean a definite 
number of systems of values of the disposable constants which will equally 
well satisfy the equations. For instance, in the theory of the second 
approximation, the equalities / 


(pr ¥ = (p2) = ($25) 
will be satisfied by supposing 7,=2,=2;*. But it is by no means evident 
a priori that this system of equalities will correspond to the absolute 
minimum of which we are in quest: nay, though even we had pa, = ¢a.=¢a3, 
those equations do not necessarily imply «,=2,=a;. Of the multiplicity of 
solutions referred to, one only gives the true minimum; but to assign à priori 
the distinguishing marks of this truest and best, hic labor, hoc opus est. It 
will be delightful to find, if it turn out to be true, that for the best form, 


g representing yX (P being a rational function of the ith degree, and Q of 
the (¢—1)th in æ), the rational quantity 


IOT WEP (Qyx) 


must be a perfect (2:—1)th power of a linear function of æ; but in the 
present state of my ignorance I dare not do more than affirm that there is a 
bare probability in favour of this being true: whoever shall first succeed in 
discovering the true form of the expression will have established a remark- 
able theorem. Here for the moment I break off, contented with having 
pointed to a theory as yet, if the expression may be allowed, sleeping in its 
cradle, but destined, I am persuaded, at no distant day to set in motion 
as large a mass of algebraical thought as has been set in motion by the 
never-to-be-forgotten Hessian discussion of the flexures of the cubic curve,— 
the turning-point between the old algebra and the new. 


Henceforward Poncelet’s theorem figures no longer as a detached method, 
a mere stroke of art in aid of the computer, but becomes integrally attached 
to the grand and progressive body of doctrine of the modern algebra. 


* If this is so, we shall have for determining the four constants the following equations : 
LEE X Ga=Gb= — r. 
But more probable than this seems the conjecture, that, supposing 2,, %a, Xg to be arranged in 
the order of their relative magnitudes, the determining equations might be 
T}=%3, pa=pb= pt, = - PX. 
Or is it possible that the character of the solution may be discontinuous, and may depend upon 
the magnitudes, relative or absolute, of the given limits a and b? Probably Dr Tchebitcheff 
would be able better than any other living analyst to answer these queries. But what an endless 
vista of future research does the prosecution of the Ponceletic method open out to us! 
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